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1. Introduction 



Observations of new physics at the Large Hadron Collider will come from the analysis of many 
scattering processes with complex final states. It is essential to prepare the theoretical groundwork 
by performing precision calculations of Standard Model production processes incorporating at least 
next-to-leading order QCD corrections. This effort requires efficient algorithms for computing multi- 
leg one-loop amplitudes. 

The unitarity method introduced in ]I] is designed to compute amplitudes by applying a uni- 
tarity cut to an amplitude on one hand, and its expansion in a basis of master integrals on the 
other H D, |], From knowledge of the basis and the general structure of the coefficients in the 
expansion, the coefficients can be constrained. 

The holomorphic anomaly reduces the problem of phase space integration to one of algebraic 
manipulation, namely evaluating residues of a complex function. By applying this operation within 
the unitarity method, coefficients can be extracted systematically. The reason this is possible is that 
the unitarity cuts of master integrals are uniquely identifiable as analytic expressions. Accordingly, 
a method was introduced to evaluate any finite four-dimensional unitarity cut and systematically 
derive compact expressions for the coefficients |7], [5[]. The evaluation was carried out in the context 
of the spinor formalism || |K], 11], [12|, |1| T4|. In [ lSj , we wrote down these general, compact 



formulas for master integral coefficients. 

The main purpose of the present paper is to improve upon those formulas in two respects. First, 



the coefficients were written as residues of the explicit formulas in |yj. Identifying the residue of 
a function simply involves performing a series expansion, but within the spinor formalism, this 
expansion is not very transparent in the case of multiple poles. Additional instructions were given 
to aid in automatizing this step. Here, our formulas will be given in terms of a truncated series 
expansion in a single scalar variable. Second, the starting point of the formulas of ]l5j , from which 
to take input data, was the result of some spinorial manipulations of the initial cut integrand. Here, 
the input data are determined directly from the initial cut integrand, as assembled from tree-level 
amplitudes. 

In this paper, we have thus eliminated the need for applying any analytic spinor identities. 
Programming the final formulas is completely straightforward. The values of the coefficients are of 
course identical to those from the formulas of our previous paper. Other general expressions for 
coefficients derived from unitarity cuts and generalized unitarity cuts of one-loop amplitudes have 
been given in 0, |17|, [18], 0, glj. 

Starting from analytic expressions for color-ordered tree amplitudes, we set up the unitarity cut 
integral. If K is the momentum in the unitarity cut, then the two cut propagators can be denoted 
by p and p — K. See Figure |]. In terms of its dependence on the loop momentum p, the cut 
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integral is a sum of terms of the following form: 



C = c I d^ 2e p 



■5 {+ \p 2 )S^\(p- Kf 



Here, c is the prefactor independent of p (but may depend on /i 2 (or -u discussed below) in our 
dimension regularization) , and the values of the momentum vectors Kj are sums of momenta of 
cyclically adjacent external particles. We work in the four- dimensional helicity scheme, so that all 
external momenta Ki are 4-dimensional and only the internal momentum p is (4 — 2e)-dimensional. 
Hence we decompose the loop momentum as |TL]. [22| 



V 



(1.2) 



where I is 4-dimensional and fl is (— 2e)-dimensional, and we further define the extra-dimensional 
parameter u by 

4/i 2 



u 

Let us then define the following four-vectors 

Qi = 

Ri = 



K 2 ' 



( /i \r , K* - {1 - y/T=u){K r K) „ 

(VI - u)K 3 + — K, 



(VI - m)P, ^ K. 




1.3) 



;i.4) 



Figure 1: Representation of the cut integral. K is the sum of external momenta on one side of the cut. 
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In terms of these momentum vectors, the cut integral may be expressed as 

c = c / *, „-~ / « *> v ^(^ Jg^gggg (1 . 5) 

where we have set n = m — k, and \£) and \£] are homogeneous spinors. This follows from the basic 
steps of spinor integration, which are reviewed in Appendix [A] but are not needed to apply our 
formulas for coefficients. The point is that the only thing we need to do is treat a general integrand 
of the form 

{VT^u){K*y +n — i npm^ (L6) 

(e\K\e] n+2 nti (t\QM\ 

The result of this integration is the subject of this paper. In terms of the vectors defined in (|1.4j ) 
from the initial data of ( |1 . 1|) , and the two integers k and n, we give formulas for the four- dimensional 
coefficients. For renormalizable theories, n < 2. Terms with n < —2 contribute to box integrals 
only; terms with n — — 1 contribute to triangle and box integrals; and terms with n > contribute 
to bubble, triangle and box integrals. To proceed to the full (i-dimensional coefficients, including 
those for pentagons, one would perform the final integral over u with the recursion and reduction 



formulas of [23, 24 



We wish to remark on a few features of our formulas. 

• Our starting point is the most general expression in field theory with a unitarity cut of a 
one-loop amplitude. Particles can be massless or massive, although in this paper we focus on 
massless propagators. Generalization to massive propagators should be straightforward. The 
propagator can be a scalar, fermion, or vector, as long as the proper degrees of freedom are 
accounted for. 1 

• The analytic formulas for tree amplitudes needed as input should be free of unphysical sin- 
gularities involving the loop momentum p, so that the form of (|1.1|) is apparent. This is 
especially important when using on-shell recursion relations to derive tree-level amplitudes. 



Using Feynman diagrams or Berends-Giele recursion [^5j to get tree-level expressions auto- 
matically circumvents this problem. 

With our formulas, we can calculate any particular coefficient directly without reference to 
other coefficients. 



^^Onc way to do this is to use Feynman diagrams to write down the full one-loop integrand expression, then 
multiply by (p 2 — m^) _1 ((p — K) 2 — m 2 ,)^ 1 along with the two delta-functions 5(p 2 — m 2 )5((p — K) 2 — m 2 ^), and use 
these to set e.g. p 2 = m\ in the integrand. In this way, even if propagators are fermions or vectors, we have counted 
everything. 
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Our formulas can easily be used to obtain the 4-dimensional part of the coefficients only, 2 by 
taking the limit u — > in ( |1.4|) . However, to be sure that all intermediate formulas will be 
well-defined, it is safest to take this limit at the end of the calculation. If we do wish to set 



it — > at the beginning beginning, some care must be taken, as discussed in Section |3T2 . 



Our formulas work for any n, although for renormalizable theories we will have n < 2. But if 
we consider (super)gravity or use a bad gauge choice, then we would have n > 2. 

In spinor notation, we will find factors of the form in the numerator. This can be 

rewritten as — 2p ■ P with P = X a Xb- So Pi can take complex values in (|1.4|). 



The formulas are presented in Section |[ In Section [3], we present an example of a 4-dimensional 
unitarity cut in one helicity configuration of the six-gluon amplitude. In Section |4|, we give examples 
from the four-gluon amplitude but keep the full ^-dimensional dependence. In Section ||, we discuss 
future applications and comparisons to other techniques. Appendix [X] reviews the first steps in 
spinor integration which form the basis of the derivation of our coefficients. Most details of the 
derivation of our formulas are presented in Appendix [B]. Appendix contains the formulas for 
triangle coefficients after the truncated series expansion has been carried out explicitly, although 
direct the use of the result in Section || is likely to be simpler. 

2. Coefficients of box, triangle and bubble integrals 

Here we give the results for the box, triangle, and bubble coefficients in the unitarity cut defined 
by the momentum K, starting from the integrand (|1 . 1| ) without the prefactor c, and using the 
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definitions ( |1.4j ). Our convention is that the n-point scalar function is defined by 

J = ^ 4vr )(4-2e)/2 f 1 , 2 ^ 

J ( 2 ^ 2e P 2 (P - K x f{p - K, - K 2 )\..(p - E;- 1 K,Y ' 

The spinor notation we use here, which may differ from other conventions, is defined as follows. 
For a four- vector ki satisfying kf = 0, 

Xi = u + (ki), \i = u„(ki), (2.2) 

thus we have the following inner products: 

(ij) = ^-\] + )=U-(k i )u + (k J ), [ij] = =u + (k i )u4k j ) (2.3) 



2 That is to say, neglecting possible rational terms which can be calculated by other methods, for example by the 
recursive techniques of |26| [30| or the specialized diagrammatic reductions of |3l]. ^] . 

3 We omit the prefactor (— that is common elsewhere in the literature O, m. 
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Note that in this paper we use "twistor" sign conventions, so that 



2k % ■ kj = (i j) [i j] (2.4) 

which differs from the standard QCD convention by a minus sign for each spinor product [i j]. Our 
definitions implying the following relations: 

(i\P\j] = u-{h) P u-fc), (i\PiPz\j) = u-(h) ft ft u + (kj) (2.5) 

The full d- dimensional amplitude will generically include pentagons in the basis. The identi- 
fication of pentagon coefficients has already been described in [24}|. The operation occurs in the 



final integral over u. Since our purpose here is to give the results of the 4-dimensional integration, 
we will not now comment any further on pentagons. In cases involving massive species, tadpole 
integrals can also arise. Unitarity methods cannot detect these. However, we expect that it will be 
possible to fix tadpole coefficients from other considerations, such as a heavy mass limit [B3 . 



2.1 Box coefficients 

A box integral is identified by the two cut propagators plus two additional ones. Following the 
setup of the previous section, denote the two additional momenta associated to a box by K r and 
K s . Then, define the vectors Q r and Q s as in ( |1.4|) . From these two vectors, we construct two null 
vectors P sr> i and P sr> 2 as follows: 4 

A sr = (2Q S ■ Q r ) 2 - AQ 2 s Ql 

p _ n , f ~2Q s -Q r + ^AT r \ 



-2Q S ■ Q r — a/As 



Sr,2 = Qs+[ ^ S Zr,2 J Qr ( 2 -6) 



Then, the box coefficient with momenta K, K r , K s is given by 



(Psr,l\K\Psr,2] U+2 Ylt=l,t=ii,j (Psr,l\Qt\Psr,^[ 



C[Q r ,Qs,K} = ( - — £■ __ . + {Psr l „ Psr2} ) . (2.7) 



4 Here we see that the formula we give is ill-defined in special cases where Q% = 0. This case can arise if we 
have set u to zero to find a 4-dimensional coefficient, and the external momentum K r is null. However, there is no 
difficulty with the underlying method. For a given box, Q r and Q s can be exchanged. Clearly, if both K r and K s 
are null, we can simply take P sr ,i = Qs,Psr,2 = Qr- ln practice, this problem can always be avoided by keeping u 
finite until the end of the calculation. 
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2.2 Triangle coefficients 



A triangle integral is identified by the two cut propagators plus one additional one. If the additional 
momentum variable is K s , then define the vector Q s as in (|1.4j) . Now construct two null vectors 
P s i and P Sj 2 as follows: 



(2Q S ■ K f - AQ 2 S K 2 



p. 



s.2 



Q s + 



2K 2 
-2Q S ■ K - 



2K 2 



K 



(2.8) 



Then, the triangle coefficient with momenta K, K s is given by 
(K 2 ) 1+n 1 1 



C[Q S ,K] 



/Al)^ 1 {n + 1)1 (P atl P a , 2 ) 



n+1 



d n+1 ( EI -=i W (Ps,i - rP s , 2 \R,Qs\P s ,i - rP St2 ) 
drn+1 ^ {PsA _ rPs)2 |g t g s | Ps)1 _ tPs2 ) 



{Ps,l <"> Ps,2} 



(2.9) 



r=0 



In practice, the multiple derivative is easy to perform in a symbolic manipulation program, ei- 
ther analytically or numerically, and we believe this is an efficient presentation of the coefficient. 
However, we have found closed expressions, and these are given in Appendix |C| for n < 2. 
If n < —2, then the coefficient is simply zero. 

2.3 Bubble coefficients 

Every unitarity cut singles out a unique bubble integral. However, in our derivation, bubble and 
triangle integrals are related, so the following formulas still require the quantities defined in ( |2.8| ). 
We further introduce two arbitrary real null vectors, 5 rj and fj, and their associated spinors. These 
null vectors must, however, be chosen generically: they should not coincide with other momentum 
variables. 

The coefficient of the bubble integral with momentum K is given by 



C[K] = (K 2 ) 1+n J2 



k n 



(r;a— q;l) , 



d<1 [ B {0) i s ) — N \ ( ft'"-"-' 
q\ ds q \ n > n ~ qy ' ^ ^ \ n > n ~ a 

q=0 \ r=l a=q 



B 



(r;a-q;2). 



, (2.10) 



'The reality condition is important. These vectors should be physical momenta of massless particles. 
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where 
1(0) 



1 (2 V -k)™ u;^ + SV )\£) 

dr n \n\[i}\T}K\q\ n {t + l)(K 2 ) t+1 (£ v ) n+1 rg =1 (£\Q P (K + srj)\£) 



bu 



, 1 (P r ,i - rP r)2 \ V \P r ^ 

/A;'"" 1 (Pr,l Pr,2) bdTb \( t+1 ) (Pr,l-rP r , 2 \K\P rtl 



7+1 



X 



(2.11) 



r=0 



(fi.i ~ rP r , 2 \Q rV \P rA - rP r , 2 ) b (fi.i - rPr,2\Rj(K + s V )\P rA - rP r , 2 ) \ 
(P r>1 - rP rj2 \r)K\P r>1 - rP r , 2 ) n+1 Y[ k p=1 ^ r (P r ,i - rP r>2 \Q p (K + s V )\P rjl - rP r , 2 ) ) t=q 1 ' 



B 



(r;6;2). 



\b+l 



t+1 



71, t 



s = 



1 (P r ,2 - rP rA \ V \Pr,2. 

t+1 



b\^K h+l (P rA P r , 2 }» dr b Ut + 1) (P r>2 - rP r>1 \K\P r;2 



(Pr,2 ~ TPrAQAP^ - rP r>1 ) b ~ tP^R^K + S7y)|P r , 2 - tP^) 

(P r ,2 - rP rjl \ V K\P r>2 - rP ril ) n+1 U h p =i,P^r ( p r,2 ~ rP r , x \Q p {K + s V )\P r , 2 - rP rjl ) 



(2.13) 



T=0 



Before ending this section, we want to make an additional remark. The derivation of these 
formulas, given in Appendix [B|, involved reducing the degree of A in both numerators and denom- 
inators. However, we could just as well choose to reduce the degree of A instead. In this case we 
would get formulas with the following replacement: |*) — > |*] and |*] — > \-k). These two sets of 
formulas are equivalent to each other in the case u ^ 0. But if we naively set u = from the 
beginning, it is possible that one of the two sets of formulas will break down. Such an example will 
be seen in Section 



3. An example from the six-gluon amplitude 

In this section we test our formulas by computing some coefficients from a one-loop partial amplitude 
with six external gluons and an adjoint scalar circulating in the loop. These contribute to the full six- 
gluon amplitude in the spinor-helicity formalism in the context of the supersymmetric decomposition 
|34] , H, |35[. The box coefficient was first computed in f36fl , and the bubble coefficient was first 
computed in || . The two- mass-triangle coefficients have not appeared in this form before, because 
it is possible to modify the basis and eliminate the corresponding integrals, as described in |7|]. 6 
In @ it was shown that for gluon amplitudes, these coefficients are constrained by IR and UV 
divergences, and we use that relation here as a consistency check. 

6 In practice it is probably best to use the modified basis, because it is less divergent. 
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Here we set the dimensional parameter u to zero from the start in order to work with simpler 
expressions. As described above, this simplification requires some care, and in fact we will see the 
consequences when we derive the triangle coefficients. 

We choose the unitarity cut of the momentum K = k 4 + k§ + fcg m the helicity configuration 
(1-2-3+4-5+6+). The cut integral is 

C 123 = J dn 2A(£r,l-,2-,3+,£+)A((-£ 2 )-,4-,5+,6+,(-4) + ) 

2 f (4 4) 2 (4 £ 2 )[3 £^[3 £ 2 \ 

" S 456[12][2 3](4 5)(5 6)7 ^ (6 4) [4 1] 

^ilC (3.1) 



S456 [1 2] [2 3] (4 5) (5 6) J n (£ - hf{£ + k x ) 
I 2{Mm (dji m ^ m ^ 



s 456 [1 2] [2 3] (4 5) (5 6) J r (£ - K)\£ + k x f 
So we have two terms, each with k = 2, and 

K x = k e , K 2 = -ki, 

so 

Qi — — Q 2 = k\. 
In the first term m = 4, and in the second term m = 3. We have 

i?i = — A1A6, R 2 = R3 = R4 = — A4A3. 

3.1 Box coefficient 

Since k = 2, we see immediately that there can be only one nonvanishing box coefficient in this cut. 
We compute the null vectors P srj ± and P sr<2 from the definitions in Q2.6Q , and define the associated 
spinors as follows. 7 

Pi2,i = Qi = -h |Pi2,i) = |6) |Pi 2 ,i] = -|6] 

^12,2 = ^2 = ^ |Pl2, 2 ) = |l) |Pl2, 2 ] = |1] 



7 We could just as well use 

Pu.i = Q2 = h, |P 21)1 ) = |1) , |P 21)1 ] = |1], 
P21.2 =Qi = -k 6 , |P 2 i, 2 ) = |6) , |P 2 i, 2 ] - -|6]. 
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Applying (|2.7| ) to the expressions under the integral signs in ( |3.1|) , we get 

;-l) n+2 s 61 (6 4) n+1 [3 lf +1 ' 



C[Q 1: Q 2 ,K] 



-2+n 
'456 



2 \ (6|K|l] n+2 
Now attach the prefactors for each of the two terms. For the first term with n = 2: 

4**61 (6 4) 3 [3 if 



S45 6 [l 2] [2 3] (4 5) (5 6) 
For the second term with n = 1: 
2(A\K\3] 



xC[Q 1} Q 2 ,K] 



xC[Q 1 ,Q 2 ,K] 



[1 2] [2 3] (4 5) (5 6) (Q\K\lf 



4^6i (6 4) 2 [3 if (4|K|3] 



s 456 [l 2] [2 3] (4 5) (5 6) — l^'«*""J (1 2][2 3] (4 5) (5 g) (p\ K \if 

So the total coefficient of the box (1|23|45|6), for the scalar contribution, is 

4**6i (6 4) 2 [3 1] 2 (6|K|3] (A\K\1] 
[1 2] [2 3] (4 5) (5 6) <6|^|if 



(3.2) 



This agrees with the expression in |36| when we incorporate the usual factor of 2/(s456S6i)- 



3.2 Triangle coefficients 

Since k = 2, we see immediately that there can be only two nonvanishing triangle coefficients in 
this cut. 



For Q 



i- 



V^l — ~" s 456 + s 45 

Pi,i = -fa 

Pl,2 



S45&—S4, 



(fa + fa) - ^fa 

S456 V * 0J S456 ° 



\Pl,2) 



|6) 

K\6] 

S45& 



\Pl,l) 

\Pl,2] 



-|6] 
K\6) 



For Q 2 



V^2 — — *456 + s 23 

p 2 ,i = fa 



S456 



{k 2 + fa 



\P2,l) 
\P2,2) 



|1) 

S456 



\P2,l] 
1^2,2] 



|1] 



-if|l) 



Let us first consider the triangle (1|23|456), with momenta K and (^2)- With the identity 



d n+1 f(a-Tb) n+1 (-r] 



dr n+1 



n+2 



c — rd 



{b — to) 
d — tc 



n+l 



(n + l)\(bc - ad) 
d n + 2 



n+l 
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we see that the formula ( |2.9|) for triangle coefficients becomes 

slt£(l\K\l} (4 6) n+1 [3 l] n+1 



C[Q 2 ,K] 



mm 



n+2 



Adding the n = 1 and n = 2 contributions and attaching the prefactors, we find that the total 
coefficient is 



*2b6<W1 ( 4 6) 3 [3 1] 



S45 6 [l 2] [2 3] (4 5) (5 6) 
2<4|JT|3] 




s i56 [l 2][2 3](4 5)(5 6) ^ 2 
s± 56 (l\K\l] (A\K\1] (6\K\3] [3 if (4 6) 2 



(3.3) 



[12][2 3](4 5)(5 6)(6|K|1] 4 

Now consider the triangle (123|45|6), with momenta K and K\ (Qi). A naive application of 
the formula ( |2.9|) in the limit u — > gives 

C[Q h K] = 

because of the R\Q\ contraction in the numerator's product. But the triangle coefficient does not 
actually vanish! This is clear, because this triangle is related by conjugation and label permutation 
to the previous one. This is the degenerate case that we discussed at the end of the previous section. 
The reason is clear. From ( p.l|) , we see that for Q 2 , the pole is [£ 1], while for Q\ it is (£ 6). Since 
the formula given in previous section is obtained by writing total derivative in [da d^\, they are 
not suitable for pole (£ 6). To deal with it we need to use the conjugate formula where we replace 
I*) — ► I*] and I*] — > I*), i.e., writing a total derivative of the form (da d\). This will be a general 
rule in all 4-dimensional calculations with null poles. We want to emphasize that this situation will 
not arise if we keep u^O until the end. 

After clarifying the subtle point we can continue our calculation. Either by taking the conjugate 
of the formula fl2.10p , which is 

'n+l / \TD ^DlD n ID ^ D 1 [ D ^DlD D ID n- TD ln+1 

{Pi 



6 123 1 


f[P 1 -rP 2 \R 1 Q 1 \ 


P 1 -rP 2 ][P 1 - 


- tP 2 \R 2 Qi\ 


Pi ~ rP 2 ] n+1 


2(n + l)!( v / A7) n+1 [A p 2\ n+1 dr ' n+1 


K [Pi-rP 2 [Q 2 Qt\ 


Pi ~ 


tP 2 ] 





or by directly applying the relabeling and conjugation to ( |3.3| ), we find 

s n +Q (6|if|6] (6 4) n+1 [1 3]™ +1 



C[Qi,K] 



2(6\K\1] 



n+2 
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Adding the two terms from n = 2 and n = 1, with prefactors, we get 



S456 (6|tf|6] (6 4) 2 [1 3f 



[1 2] [2 3] (4 5) (5 6) 
g456 (6\K\6] (6 4) 2 [1 3] 2 (4\K\1] (6\K\3] 
[1 2] [2 3] (4 5) (5 6) (6|^|1] 4 



s 456 (6 4) [1 3] + <4|*T|3] (6|#|1]) 



(3.4) 



Consistency check: 



These particular triangle coefficients have not been isolated before, because two-mass triangles 
disappear in the modified integral basis proposed in J7[. We can now perform a consistency check 
based on the same identity that allowed the basis to be modified. Consider all the contributions 
to the divergence (— s)~ e , where s = K 2 . In these example, there are exactly these two 2-mass 
triangles plus the single box from the previous subsection. The condition expressing the vanishing 
of this divergence is 8 

2 V V / . x 1 







J2m h 



St 



-t) 



(3.5) 



2m h 



C[l|23|456]- 



C[123|45]6] " 



2m h 
C 4 



<S456 S 61 ' ' ~ s 456 + s 23 

2 1 1 

J^sV, + C[1|23|456] (iWJ ~ C[1|23|456, (W| 



S456 + s 45 



It is easy to see that this identity is satisfied by our coefficients given in fl3.2|) ,(|3T 



3.3 Bubble coefficient 



4. This choice gives somewhat simpler formulas; for example, 
B y "l_ n {s) and B^'"Z^ 1 '' '(s) are identically zero, because there is a sufficiently high power of r inside 



Let us choose rj = 3 and fj 
^ n _ q (s) andi^-^t 
the derivative. We find 



B 



s) 



:-!)« [6 3] 



a—q 



d a ~ q 



[3 6] 



n—a+l 



(6\K\6] 



n-q+2 



q)\dr a - 



n+l-q, 
5 456 1 



n 



a + 1) 



x 



r a -«(s 456 [6\K + sk 3 \6) - rs [6\3\K\6])( S4BB (4 6) - r (4\K\6\) 



n+l 



(s 456 (3 6) - r (3\K\6})i(s A56 [1\K + sk 3 \6) 



T 



\K + sk 3 \K 



T=0 



^The relative sign between box and triangle terms comes because our sign conventions for the master integrals 



(2.1) differ from those of 



12 



(l;a~q;2) f 



;-!)« [6 3] 



a—q 



d a ~ q 



[3|A-|6> 



n— a+1 



|K|6] n - g+i (a - g)! dr«-<? ^^"'(ti - a + 1] 

(s [Q\3\K\Q] - T5456 [6\K + sfc 3 |6))((4|A:|6] - rs 456 (4 6))" +1 
((3\K\6] - r^e (3 6>H[1|K + sfcsl^le] - r S456 [1|AT + sAs 3 |6» 



T = 



B%a ;1 \s) = o 



»(2;o-g;2)/ \ 

*"*n,n— a v" 5 / 



-1)" [1 3] a " 



[3|tf|l> 



n—a+l 



{l\K\l 



in— 3+1 



q)\ dr a - 



n+l-g 
3 456 



(n 



((4|ir|l]-TS456 (4 1))' 



n+1 



((3|tf|l] - rs 456 (3 l))*«6|tf|l] - ts 456 (6 1)) 



T = 



We then substitute these expressions into (|2.10f) and attach the prefactors. For (s), in fact 

only the q = contributions matter, because the s-dependence has dropped out with our choice of 
\rj) = |3). We have checked numerically that the result agrees with the corresponding result derived 
by the technique of M. 9 



4. A d-dimensional example: four gluons 



In this section we illustrate the use of the formulas in Section ^ in the case of four gluons with 
a scalar propagating in the loop. These amplitudes were first given in p2| . Our notation and 
presentation here are more similar to J57J] and especially pi ], where these amplitudes were derived 
by newer techniques. Here we have verified that our results reproduce those in the literature. 10 
We stop just before the final integral over u, which could in general be done by the techniques of 
P3j UM , developed in the context of rf-dimensional unitarity ^2], [J7|. Note therefore that 

the labels of "box, triangle, bubble" are used in the (i-dimensional sense. The four-gluon amplitude 
is a nice test of our formulas, because this simple case is where they are most likely to break down, 
for example by a bad choice of rj, as we shall see in the last configuration. We consider three of the 
four independent helicity configurations, since the fourth adds no new features. 
For ease of presentation, we make use of the variable z as given in (/O). 



9 We were unable to confirm numerically the printed value of the corresponding coefficient in so we repeated 
the calculation. 



10 



The — I h amplitude was given the wrong overall sign in equation (4.31) of 



13 



4.1 (l+,2+,3+ 4+) 

The simplest helicity configuration is (+ + ++)• The integrand for the cut K = K\ 2 is 

2//[12][3 4] 1 
(1 2) (3 4) (p-ktfip + k*) 2 

From this, by comparing with general formula ( |l.lj ) we have m — 0, k = 2, K\ = — — ^4, 

thus n = m — k = —2, so there are neither triangle nor bubble contributions. There is only one box 
coefficient. There are no vectors Pj. The expression inside the parentheses in fl2.7|) thus degenerates 
to 1, so we find 



2 / /[12][3 4] l 2//[12][3 4] 

(1 2) (3 4) 2 l ' (1 2) (3 4) ' 1 ' ; 



(1 2) (3 4) 2 V y (1 2) (3 4) 
4.2 (l-,2+,3 + ,4 + ) 

For this case we have two cuts C\ 2 and C 4 i. These two cuts are related to each other by symmetry, 
so we focus on cut C41. The integrand is 

u[2 3] <lb|4] 2 



2(2 3) (p-k 4 )(p+k 3 y 
For this case we have m = k = 2 so n = m — k = 0, Ki = fc 4 , = — k 3 , thus 

P = A!A 4 , P = i?! = P 2 = -(1 - 2z)AiA 4 
g 1 = - (1 - z)k 4} Q 2 = (1 - z)/^ + zfc 2 



Box: Using 



Q 2 = Ql = -K\ x , 2Q X ■ Q 2 = -K\ 2 + |(# 2 2 - ^ 2 3 ) 
A12 = s 2 (l -u)(l + u-), s = K\ 2 , t = fT 2 3 , 



we have 



u[2 3] [4|3|1) 2 (2 + An)s 2 1 s 41 [2 3] 2 [4 3] 2 w(2 + Au) ^ = £13 (4 2) 



2(2 3) 2s 2 3 2 8[1 3] 2 ' s 12 

Triangle: 

Formally, there are two triangles identified by the momenta K\ and K 2 . In this special example 
with only four external particles, they are, in fact, the same triangle. So we will need to add these 
two contributions to the final coefficient. 
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Let us start with K 1 . From Q2.8Q , we find 

A 1 = (l-2z) 2 4, P ltl = {l-2z)k 1 , Pi, 2 = -(l-2z)k 4 . 
For the spinors, we choose 

\Pi,i) = 11), lA,i] = (1 - 2^)11], |P li2 ) = |4), |P li2 ] = -(1-2*)|4]. 
Since n = 0, using Q2.9|) we get 



(1 - 2z) d / z 2 (4 l) 2 [4 l] 2 (1 4) 2 zV (4 l) 2 [4 l] 2 (1 4) 2 \ s 12 s 2 Al 



2(41) dr \ih-Th\Q2Qi\h-Th) ' (h - Th\Q 2 Qi\h - rh) J 2(4|3|1] 2 
With the prefactor included, we have 

M[23] -qg^ 41 ]- US4l[23]2[34]2 



2(2 3) ^ 4iJ 4s 12 [3 1] 2 
For the triangle with K 2 we have 

A 2 = (1 - 2z) 2 s 2 23 , P 2A = -(1 - 2z)£; 2 , P %2 = (1 - 2z)fc 3 . 

Similar calculations give 

«[2 3] ^ ^ , _ u[2 3] 2 [4 3] 2 /, 2 3 



-C[Q 2 , K A1 \ = J ro l 119 J - 2s l3 -s 12 ). 



2(2 3) ^ 4 S23 [3 1]2 V^i2 

Adding these two contributions together, we find that the triangle coefficient is 

u[2 3] 2 [4 3] 2 s 12 (g 41 - gig) 
2s 23 s 12 [3 l] 2 

Bubble: The formula ( |2.10| ) reduces to 



(4.3) 



rjr^i _ mm 2 j ^ 2 (gsiMgd (fi,il^l^ 2 ] 2 f i 

[ J [nlQiKlvMQ^] ^ l \^ r {P rA \K\P rjl }{P r MPr,2]{Pr,l\Q2\Pr,2] W ^ 

Taking 7] = it is easy to see that both the first term and the r = 1 terms are zero, so we are left 
with the r = 2 terms only. The result is 

[4 3] 2 si 3 / s 13 \ 
[3 1] 2 s 4 iV W" 
With the prefactor included, we find that the bubble coefficient is 

«[2 \ C[R] = us 13 (s 12 - s 13 )[2 3] 2 [3 4] 2 ^ 



2(2 3) 1 J 2s 2 lSl2 [13] 2 
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4.3 (l-,2+,3-,4 H 



The integrand for the cut K41 is given by 

2(1|£|4] 2 (3|£|2f 



s\M-h) 2 -^m + h) 2 -^) 

For this case we have m = 4, k = 2, so n = m — k = 2, K\ = fc 4 , K 2 = — k 3 , thus 

Ri = R 2 = -(1 - 22)AiA 4 , R 3 = R 4 = -(1 - 2z)X 3 X 2 . 
Qi = —zki - (1 - z)k±, Q 2 = (1 — + -2^2- 

Box: By now it is straightforward to find that the coefficient is 



2 rrn n ^ (13^s|i(8£|2 + 8s 12 s 13 u + s 2 13 u 2 

—C[Qi,Q 2 ,K\ = — 

s 4 i 8 ( 

Triangle: For the first triangle, with Ki, we have 

1 c? 3 / (fci — Tk^\R\Q\\ki ... ' ' ■ \ \ '■ ' -' : >' c 1. 1 ' ' 1. ' " 1 / 1, 
G^ijAj = — — , . 7^ jz , „ 1, ^ l fc 



1 d 3 


' (ki — r&; 4 


R1Q1 


fci — rA; 4 ) 2 (^i 


— r&; 4 


R3Q1 


k\ — r/c 4 ) 2 


12(l-2z) 3 (l 4) 3 rfr 3 




(A 


1 - t£;4|<22<2i 


ki — rk^) 



In this case, the factor [k\ — rk^\R-\Qi\ki — rk^) is proportional to r 2 , so the contribution of first 
term is zero and we have 

4- C [Q 1 , Jf | = -< 13 > 2 ?' 1 + I » 2 + I "), 1=£H 
si, K ' 1 2(2 4> 2 A 3 S12 

By symmetry, for the triangle with momentum K 2 we just need to do exchange the labels 1 «-» 
2,3^4, so the final coefficient is just twice what is written above, namely 

(13) 2 s 13 (l + l) 2 j2 + l M ) ~ =£l 3 , , 

(2 4) 2 A 3 ' s 12 ' 1 ' ] 

Bubble: To present this example analytically, we choose 77 = /c 4 rather than a more generic 
value. However, we do run into a problem then, because of an accidental degeneracy of poles. This 
problem could have been avoided by choosing a generic 77, but for analytic purposes, we considered 
this case separately, and it is presented in Appendix p. 3.1 , The formula ( B.28| ) for the bubble 
coefficient now takes modified input, as given in (|B.25|) , ( B.26|) and ( B.27|) . 

First, consider the term ( |B.25| ). Since n = 2, there is a third derivative in r. Since the factor 
[4 1 (if — Trf)Ri(K + ski)(K — rrf)\A\ is proportional to r 2 , this contribution vanishes. Similarly, 
( B.26|) will be zero. The only remaining part is the r = 2 case, ( B.27|) , with 



P 1 = -(1 - 2z)k 2 , P 2 = (1 - 2z)K 3 , A r=2 = (1 - 2z)s 



2 

41 • 
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Let us discuss B^f' 2 \s) first. Notice that 

(3 - t2\R 3 (K + S 7])\3 - t2) = (1 - 2z)t (2 3) s 12 B(-s + =(1 + (1 + s)A)). 

B 

We see that to get a nonzero value of B^l (s) we must have a = 2, and more specifically, 

£§ a;2) ( S ) = 0, a = 0,1 
,(2;«=2ii),„ , _ (1 - 2z)[i l?B 4 (1 2) 4 (- Sl2 )* +1 *V 



s!i(2 3) 2 (t + l) S *+ 1 (l-2z)-^) 

Because the factor s 2 appears in B^t* - ( s )> there is a nonzero contribution only when we take 
g = 2 in the derivative with respect to s. But then a — q = 0, and since B 2 2 ' a 9 ' 2 ^(s) is only nonzero 
for a — q = 2, the contribution from this part is zero. 
For we have 

(2;a;1) (1 3) 2 (-l) a (l - 2zf- a A t+l d a ( (1 - r) m (l - z + zAtY(sAt - 1 - (1 + s)Af 

B21 ' ( s ) 



(2 4) 2 A 2 s 4 i a\(t + Vjil + A)^ 1 dr a y (1 - r) 4 -«((l - 2z) - zs) 

where r = r/S, B = (2 4) / (3 4), and 1 = si 3 /s 12 . Summing up (t, a) = (2, 0), (1, 1), (0, 2) with 
s = 0, (t, a) = (1, 0), (0, 1) with the first derivative of s, and (t, a) = (0, 0) with the second derivative 
of s, we finally find that the coefficient is 

2 cm _ 2 g 3 ) 2 g i3 (12 + 3l(6 + rQ + 1 2 (4 + 5m)) I= £i3 f47 , 
si. (2 4) 2 s 4 i 12A 2 ' S12' 
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5. Discussion 



Since the formalism described here is based on unitarity cuts of the amplitude, it shares with other 
unitarity-based approaches 11 the property that the input required is simply a collection of tree-level 
amplitudes. These are manifestly gauge invariant and can take quite compact forms. By dealing 
with different cuts separately, we can attack the problem in stages. 

Furthermore, our formulas separate and identify the coefficients of individual master integrals. 
A single unitarity cut yields, directly and separately, the coefficients of all the master integrals with 
the same cut propagators. Any single coefficient can be targeted individually, without the need to 
first compute any others or additional spurious terms. 

11 For a review, see Section 4 of Ell. 
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5.1 Comparison with other approaches 

The reduction algorithm of Ossola, Papadopoulos and Pittau (OPP) J42], |3|, [44| produces coefficients 
through algebraic operations at the integrand level, through recursive solution of a set of algebraic 
equations. In fact, our formulas given here are the results of solving algebraic equations in a 
different style. In the OPP method, several points of phase space are used, while in our method, 
we differentiate at a single point of phase space. The derivative operator can be interpreted as an 
algebraic procedure as applied to rational functions at a single point. 

Just like our result, coefficients from OPP method can be fed into the <i-dimensional unitarity 
program as described in [F2^, ^M. Alternatively, the algorithm may be interpreted numerically, and 
in fact such an implementation has now been given jEJ (see also the procedure of [SJ). We believe 
that the formulas of the present paper are also well suited for numerical programming, and this will 
be the subject of forthcoming work. 

One final note on comparison to the OPP method is that our formulas are valid for arbitrary 
values of n, in particular for n > 2. Such an extension has been mentioned within the OPP method, 
although details have not been worked out. 



An approach that is closer in spirit to ours was given by Forde in [19[]. There, coefficients for 
boxes, triangles and bubbles are given within the spinor formalism. The foundation there consists 
of generalized unitarity cuts, namely quadruple cuts for boxes, triple cuts for triangles, and ordinary 
double cuts for bubbles. The motivation was to capitalize on the efficiency of quadruple cuts for 
box coefficients, and also to be able to target specific coefficients. Forde's final formulas resemble 
ours in that they are based on data from tree amplitudes and given in terms of a coefficient in 
a series expansion of one variable for triangles, and two variables for bubbles. The formula for a 
bubble coefficient, however, depends on tree amplitude input for all possible triple cuts, while ours 
comes directly from the ordinary double cut (though still depending on all possible momenta from 
a hypothetical third cut). If the aim is to assemble an amplitude in its entirety, then of course the 
complete tree-level amplitude input will be available anyway. 

In WB, we discussed the application of quadruple cuts to box and pentagon integrals in d 



dimensions. A <i-dimensional analysis of triple cuts for triangle integrals has been given in . Let 
us now briefly examine the triple cut in the context of the present paper in order to make contact 



with the result of [119 



With three delta functions in d dimensions, we first use two of them to set up the four- 
dimensional spinor integrand, as explained in Appendix |A|. After integrating over the variable 
t, defined in ( |A.2| ), we arrive at the integral 



(£ di) [£ d£] G{£) 5 



K 2 



(£\K\i\ 
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Now we can use momenta Q, K to construct two null momenta as in ( |2.8|) and expand our spinor 
variables in the basis of their spinor components as follows: 

\i) = \p x ) + z\p 2 ), \e} = \p 1 } + z\p 2 ] 

Here z is a complex number and z is its conjugate. With this substitution, we get 

dzdz(2P P)C(z~)5(k APMPi] + Z * {P21QIP2] } (5 1) 

dz dz (2P, P 2 ) G(z, z) 5 \K {miPi]+znP2lQlP2] ) ■ (5-1) 

Now we can change to polar coordinates so that z = re ld and dzdz = rdrdO. Furthermore, if we 
now define the new variable t = e l9 , we have 



dzdz = rdr x — ■ — . 

t 

The delta function depends only on r, so we can use it to integrate over r. Then we are left with 
t integration only. From here, for example, it is easy to see the vanishing condition given in eq. 
(4.20) of ||. 



Furthermore, for box integrals, we have an extra propagator, so the general form of the integrand 
is 1 / (a + tb + t -1 c). Only polynomial terms correspond to the triangle contribution. 

The parametrization we have used here is not exactly the one used by Forde, but the central 
idea is the same and t is the angle variable for both triple cuts and double cuts. 

5.2 Prospects 

The most obvious and immediate application of our results will be to the computation of complete 
one-loop amplitudes, as in the example of Section [|, where the w-dependent expressions for coeffi- 
cients are fed into the reduction formulas of []23], 24j to give the final e-dependent coefficients. The 



four-momenta of the external particles may be numerical at every step. The reduction formulas 
currently require analytic expressions in u. 

Our formulas may also be specialized to the cut-constructible part of the amplitude, as in the 
example of Section [|, simply by setting u — > at the end of the calculation and interpreting the 
formulas as exact coefficients of 4-dimensional master integrals. 

Since several methods proposed in the literature are specialized for computing either cut- 
constructible or rational components of an amplitude, it would also be very interesting to specialize 
our formulas to isolate the rational part of one-loop amplitudes. This could be done by studying 
the e-dependence of the reduction formulas together with the u-dependence in the coefficients, in 
order to focus precisely on the e° term in the final e-expansion of the amplitude. We will return to 
this point in a future publication. 
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Finally, as we remarked in the introduction, our formulas apply to amplitudes with massive or 
massless propagators. In order to arrive at complete amplitudes in the massive case, the master 
integrals should be evaluated explicitly, 12 and our results will need to be supplemented with the 
contributions of tadpole and massless bubble integrals. 
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A. Setting up the cut integral 

In this appendix we briefly review the first steps in spinor integration, in the context of <i-dimensional 
unitarity, leading from equation fll.lp to equation (|1.5|). For a fuller discussion of this technique, 
24]. Within the four- dimensional helicity scheme, we apply (|1 . | ) and ( |1.3|) . In the integrand, 



sec 



p is replaced by £, and the measure is transformed as follows 

d 4 ~ 2t p f dH (4tt)' 



dp (// 



-l-e 



(2tt) 4 - 2£ J (27r) 4 r(-e) 

We now drop the factor ^ 2 -l)^T(-e) (t") ' w ^i cn ^ s universal and common to cuts of amplitudes 
and master integrals. Following |2B], we further decompose the 4-dimensional momentum into 
a null component and a component proportional to the cut momentum K. 



While changing the variable p to u with (|1.3Q , we note that the kinematics of the unitarity cut 
constrain the integration domain to be u G [0,1]. Our cut integral (|1 . 1| ) can now been rewritten as 
the following expression. 

J du u^ e J dz (1 - 2z)S(z(l - z) - |) 
J dH 5 + (f)5((l - 2z)K 2 -2£-K) : ^=^~ 2Pi ' ^ + zK V 



Y\"{K]-z{2K r K)-2t.K j 



12 For a uniform mass, this has been done in fE2| 
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Notice that the z-integral can now be done with the first delta function. In fact, the kinematics of 
the unitarity cut require us to choose exactly one solution for z. If we take K > 0, then 

1 — a/1 — u , 

z— , or equivalently, 1 — 2^ = yl — u. (A.l) 

Now we change to spinor variables with |47|] 

t = tXX, (A.2) 

where t takes nonnegative real values, and A and A are homogeneous spinors. The measure trans- 
forms as 

dH 5 {+ \t 2 ) (•) = / dtt (A dX) [X d\](m). 

Jo Jx=x 

The domain of integration of t is again consistent with the kinematic region of the unitarity cut. 
From here on, we use \£) and \£] interchangeably with A and A. We have now arrived at the following 
expression: 

c = / *, / { e M) V wjttw- + 1 

Finally, we use the remaining delta function to perform the integral over the variable t. With the 
substitution ( |A.1[ ), the result is equation (|1.5|). 



B. Derivation of the formulas for coefficients 

In this appendix we outline the derivation of the main results of this paper, which are the formulas 
Q2.7|), ( |2.9| ) and ( [2.10| ). Our technique is the type of spinor integration carried out in [|7|, p|, p!5 



but we stress that understanding these techniques is unnecessary for applying the results. Indeed, 
equivalent formulas have already appeared in [T5[]. 13 The difference is that our starting point is now 
the raw unitarity cut integral, before converting the loop momentum to spinor variables. In the 
final formulas, we now explicitly evaluate the residue at multiple poles. Additionally, the present 
versions of the formulas feature substantial simplification of the bubble coefficients. 

Our foundation here is the framework laid out in [|15 and its references. Let us briefly recall 



the key ideas. The general integrand given as the starting point in |L5 



is 14 



term <^r +2 nLw^]' (R1) 



13 



The formulas for coefficients given in ]f 5| differ from the ones given here by a factor of yl — u. This comes from 



the convention of our starting point (1.1) or equivalently fll-q ), where this factor appears explicitly. 
14 We have redefined the index n for consistency. 
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Comparing with the expression ( |1.5| ), we see that we will take G(X) to be constant and [aj\ = (£\Rj\- 
The idea of spinor integration is to rewrite the integral so that we can carry it out with the residue 
theorem. The next step, therefore, is to isolate poles by splitting the denominator factors with 
spinor identities such as 

M _ [a\Qi\£) | NW (R2) 



(£\Qi\A(AQ2\i] (£\Q2Qi\£) (i\Qi\i] (AQiQN) (AQM 

This procedure is applied to the amplitude on one hand and the master integrals on the other. By 
matching functional forms, we extract the coefficients. 

B.l Box 



The formula ( |2.7D for a box coefficient is trivially related to the one given in ||15|| . We only need to 



observe that now that we take G(X) to be constant and [aj\ = (£\Rj\ in ( |B.1| ), the poles from the 
factors (£\Q s Q r \£) are inserted into [dj\ as well when we evaluate the residue. 

B.2 Triangle 



For a triangle associated to momenta K S ,K, the coefficient was found in [nj to be the difference 
of the two residues from the poles in (£\Q s K\£) n+2 of the following function: 15 

(-i)"(^ 2 ) 1+ "VAl Y^mQM) 

2 (£\QsK\i) n+2 U.Ut^(AQtQs\£) ' 

The quantities A s ,P Si i,P Sj2 were defined in ( |2~8|) specifically to deal with the factor (£\Q S K\£) 
by identifying the poles explicitly. With those definitions, we find an identity that separates the 
two poles: 

(£\Q S K\£) = (£ P sA ) (£ P s , 2 ) R2[Ps ± Ps ' 2] . (B.4) 



Now consider the residue from a multiple pole, in an expression of the form 

1 N(\£),\£}) 
(£ V) n D{\£) ,{£})' 

We can start by substituting \£] — \r)], so we are dealing with the holomorphic function 

i N(\£)Av}) 

(£ V ) n D(\£),\v}Y 



15 



Here we have again redefined n and made substitutions for G(A) and [a.j 
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For an arbitrary auxiliary spinor (, we have the following identity. 



d 



n—l 



(£ ( V - rC)Y 



dr n - 



(n -!)!(£ V (v-rQ) 



(B.5) 



Thus we find 



1 



d 



n—l 



1 



1 



n(\£),\ v )) 



(£ ( V - rC)) n dr^ \( n - 1)! (£ 0^ (* (V ~ <)> D(\£) , \rj\) 
Now we extract the residue at the single pole (£ (77 — tQ) before taking the derivative. We find 

d n ' x ( 1 N(\ V -t(),\v]) 



(B.6) 



dr 



in 



(B.7) 



To obtain the residues from the factor (£\QK\£) ri , we use equation ( p.4j) to rewrite it in terms 
of two multiple poles. 16 Then we apply ( |B.7| ) to compute the two residues as follows: 



Ri 
R2 



d 



n-1 



dr™- 1 
d n ~ x 



1 



[n - 1)! (Pi Ci) n_1 (^Y\Pi P 3 ] n (Pi - nCi, P 2 ) n D(\p 1 - nCi) , IAD 

1 (K 2 ) n N(\P 2 -T 2 ( 2 ),\P 2 }) 

(n-l)! (P 2 C2} n_1 (VS)*[Pi P 2 ] B (A - T2C2, A)" D(|P 2 " r 2 C 2 ) , \P 3 ]) 



where we can choose different auxiliary spinors Ci>C 2 f° r the two poles. To simplify further, we 
choose |Ci) = \P 2 ) and |( 2 ) = |Pi) and use the identity [Pi P 2 ] (Pi P 2 ) = -A/K 2 . Finally we find 



Pi 
P 2 



(-1)" 1 d n ~ l f N(\P 1 -r 1 P 2 ),\P 1 }) 

yS)» (n - 1)! (Pi P 2 ) n ~ l drr 1 \D{\I\ - r x P 2 ) , |P]) 
(_!)» 



1 d™- 1 


^iV(|P 2 -r 2 Pi), 


ft]) 


l)!(PiP 2 ) n - 1 rfr 2 "- 1 


VP(|P 2 -r 2 Pi), 


ft]) 



(B.8) 
(B.9) 



Using ( B.8|) and (|B.9|) with our original expression ( B.3|) , we get the formula (|2.9|) for the triangle 
coefficient. That formula may look as though it is not completely explicit because we still need to 
perform a differentiation. But this is easily done in a symbolic manipulation program. We do offer 
explicit formulas in Appendix |C| but do not expect those to be more useful. 

Recall that the final result must be a rational function, so the square roots from a/A7 should 
eventually combine into polynomial expressions. 
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Throughout the rest of this derivation, we drop the subscript s to avoid cluttering the formulas. 
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B.3 Bubble 



Our derivation here parallels the one in |15|], but the splitting identities are more systematic and 



the final formula is now written explicitly. 

First, we would like to split the denominator factors in ( |B.1[ ) using the following generalization 
of (BJ): 

L/mirf (R10) 



This formula is applicable when and only when all Qi and K are different. To use it, we deform 
( |B.1|) by introducing small independent parameters Sj, i — 1, n + 1 and a real null vector r/. 



G(A)n;il% 



m + SlV \£]uU {t\Q P \£\ 



(B.11) 



The final result will be recovered by taking the limit Sj — > 0. 

Now we can apply ( |B.10| ) to ( [B.ll| ) to find the following expression: 

1 G(X)U r ;= k 1 [a J \K + s lV \i) 

We can see that the Sj — > limit is smooth in the second line, resulting in terms of the form 
Fi(X) / ((£\K\£] (£\Qi\£]). We know from f| that these terms yield pure logarithms, in this case for 
the triangle integrals associated with momenta K and K^. So we restrict our attention to the first 
line, ( p,12|) . Rewrite it as 

1 G(\)Yi;+ k A^\K+s^\l) 

^ (£\K\£] (£\K + Si r,\£] (£\K v \£) n U^d 8 * - s q ) UU (Wp(K + s iV )\£) 



Now we can must take the limit Sj — > carefully. We find that the bubble coefficient is 

(-l)grfg£ W|W _ g ( a ) 



(is-? 

where we have defined the function 



(B.14) 



s=0 
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The fact that ( B.14 ) represents the bubble coefficient can be proved by induction. The case n = 
is trivial. Assume that it is true for n, and let us now introduce a single parameter s to rewrite 
PUD as 

g(\) n?-f +1 K- § 

(WT +2 m + sr ] \£]U k p= AW] 

We now treat the factor (£\K + 7rj\£) on the same footing as and apply the result for n. 

The bubble contribution can be expressed as a sum of two terms ii and 1%. The first term is 

i G(A)nS +1 K-i^+^> 



\K\£] (£\K + s V \£} (£\K(K + s V )\£) n+1 nj=i (t\Q P (K + srj)\i) 



* n+1 v t <^i^ + V (^i^K} n+1 np=i(^iQp(^+^)K) 

After taking the 7 — > limit, we have 



n+1 n+l-a , is n+1+t , a . 



Jl Z^ Z^ 5h+i-t-« a! ds° ' 1 ' 

a=0 t=0 



where B ntt (s) is defined by (|B.15| ). 
The second contribution is 



(-ljggg^g = 0) 
g=0 * 



where 

tin AS) = 



£\y\£f G(\)U;^ +1 [a,\K + s V \£) 

K\£] 2+t (£\ v K\£) n (£\ (K + 7 V ) (K + s V ) \£) UU Wvi K + s v) l*> 

' -B n+1>t {s). (B.18) 



[s — s / 

We must take the s — > limit before s — > 0, so first we substitute 

-6 



d«B nin „ q (s) 



ds<* 

to find 



gi dP- 



(-l)g rf^ 
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or equivalently, 

n n ~ a (_-\\n— t ja 

^ = EE^^W-0) (B.19) 



a=0 t=0 



Now it is easy to see that Ii + I2 is nonzero only if a + t — n + 1. Therefore we can write 

h + h = J2 ^^^Bnwx-ais = 0), (B.20) 



a\ ds a ' 

a=0 



and thus we have proved the formula ( |B.14|) for n + 1. 

Now that we have established that the bubble coefficient comes from (|B.14|) with the definition 
( |B.15p , we need to identify the poles and find the residues. 

Rewrite the integrand ( p,14|) as a total derivative by using 

Now let us specialize to the integrand of ( |1.5|) , so that G(X) is constant and [aj\ = (£\Rj\. Now we 
define 17 

b. M s _JL_ <MT rSgmfeW 

A (« + l)(<|A'K]' +1 (<'|r,A'K)" +1 n*,i(« P (A' + S »))l<'> 

Here it is important that 77 be completely generic, so that there are no accidental degeneracies. For 
an alternative approach, see Subsection |B.3.1| . There are three kinds of poles. The first, at i = 77, 
has no residue because the numerator factor (^M^]*" 1 " 1 becomes zero. The second, at \£) = K\rj\, is 
a multiple pole of the type discussed in [B.2| , so we see that its residue is ( |2.11| ). The last kind of 
pole is from the factor (£\Q r (K + sn)\£). Here we perform a series expansion in the parameter s, 
which we will ultimately set to zero. The expansion is 

(£\Q r (K + s V )\£) ^ ] (£\Q r K\£) a+1 ' 1 ' ] 

The residue is then B ^f' (s) — B^[f' 2 \s), with the definitions given in (|2.12|) and Q2.13Q . 
Combining these contributions, we find that the sum of the residues at poles of B n> t(s) is 

sSoo + EE(- s ) a - b^\s)) . (B.23) 

r=l a=0 

Feeding (P-23| ) into ( |B.14|) and simplifying the result gives us our final expression for the bubble 
coefficient, (|2T0D . 



17 



The B n t (s) is the splitting result while /3„ it (s) is after writing into total derivative, i.e.,B n t (s) = [dX d-^\B n ^{s). 
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B.3.1 A special choice of r] 

In this appendix, we describe the consequences of choosing rj = K\ in the case where K\ . This 
choice may be convenient for small examples worked by hand, but we emphatically recommend 
choosing a generic rj wherever possible. 

The reason that such a special choice of i] presents a problem is the following. From ( |1.4|) , we 
can see that 



\QiK\ 



-(1 - 2z) (eiKiKlQ = -(1 - 2z) (£\r]K\£) 



Therefore, in the expression ( |5~2lT ), the poles from (l\r]K\£) and (£\Qi(K + sr))\£), will overlap, and 
the way we read off their residues should be modified respectively. 
In this special case, it is easy to see that 



1 



{£\Q X {K + s V )\£) 
Now instead of (|B.21| ), we have 

1 



£<-< 

a=0 



J£\QMTf_ 
(t\QiK\£) a+l 



1 



1 



(£\VK\£) (1 - 2z) 



sz 



K 2 



B n , t (s) 



(£\ V \£] 



t+i 



W2 + s V )\£) 



(1 - 2z) - sz&P (t + 1) (£\K\£] t+1 (£\ v K\£) n+2 U k p=2 (i\Q P (K + s V ) \£) ' 
Continuing this way, we find 

d n+1 ( 1 [ritfKlr]]-^ 1 f (2,/ ■ /\ j x 



dT n+i l(i _ 2z) — sz^fi- (t + l)(n + 1)! 



K 2 

(* V) n+2 UU (£\Q P (K + srj)\£) 



K 2 



\i)^\K-rrj\ V ] 



(B.24) 



(B.25) 



B 



(r;a;l) 



n.t 



{-If 



d a 



t+i 



/ ST +1 a! 



x 



1 - 2z) - sz^j 

y\£) a U-=i^\ ^iK + srm \ 
(£\Q P (K + sr})\£) J 



(Pr,l Pr,2) a ^ \( t + 1 ) (£\K\ 



it+1 



{£\vK\£) n+ 'X[ 



\q = \Pr,l],\t) = \Pr,l)-T\P r ,2) 



[B.26) 



B 



(r;a;2) 
n,t 1 



1 



d a 



s = 



t+i 



(1 - 2z) 



sz 



K 2 



^A~ r a+1 a\ (P r>1 P r , 2 ) a dr a Ut + 1) (£\K\ 



it+i 



x 



(t\Qrv\£) a n;^ + s V )\£) \ 

{£\rjK\£) n+2 UUv^r (i\Q P (K + s V )\£)J 



\i] = \Pr,2]M) = \ P r,2)-r\P r ,l) 



[B.27) 
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and the coefficient is given by 



c[ K ] n = (KY +n E (b^s) + e E (B { ::::f\s) - <rr 2) (*; 

g=0 ^' S \ r=2 a=q 



. (B.28) 



s=0 



where the definitions of the functions B must be taken from ( |B.25|) , (|B.26| ) and (|B.27|) . 



C. Closed forms for triangle coefficients 

We have given the general expression for coefficients of triangles as a formula involving a multiple 
derivative (|2.9| ). We can also carry out the differentiation explicitly. For practical purposes, we 
need to consider only cases with n < 2. 

When n < —2, the contribution is simply zero. 

When n = — 1, there is no derivative, so the result is just 



C[Q s ,K} n -_ 



2 vnt Ws (^ii^i^2 

When n = it is given by 

k 2 \ nli^ii^i^] 



(C.l) 



C[Q 8 ,K] n=0 



2A - \nlu*,(p.,i\Qt\p.,2] 



E 



(2Q s -K)(2R r Q s )-2Q 2 s (2R r K) 



- E 

t=l,t^s 



{2Q s -K){2Q t -Q s )-2Q 2 s {2Q t -K) 



\Ps,l\Qt\Ps,2 

When n = 1 it will be 

{K^f r nS(^,ii^i^, 2 ] 



\-Ps,l|-Rj|-Ps,2] 



{Ps,l <"> Ps,2} 



(C.2) 



4A " {Ut^s(Ps,l\Qt\P 



'fc+1 

E 



(2Q 8 -K)(2R j -Q 8 )-2Q 2 8 (2R r K) 

(P.,l\Rj\P.,2] 



(2Q s -K)(2Q t -Q s )-2Q 2 (2Q t -K ) 

(P S ,l\Qt\P S ,2\ 

[(2Q S ■ K)(2R j ■ Q s ) - 2Q 2 S {2R 3 ■ K)} 2 + 2Q 2 if 2 (P^i^P^] (P^lP^] 



- E 

fc+1 

+E : 

j'=i 
fc 

- E 

t=l,t^s 

+{P SI ! <-> p s , 2 }} 



-[(2Q S ■ K)(2Q t ■ Q s ) - 2Q 2 s (2Q t ■ K)} 2 + 2Q 2 S K 2 (P s ,i|Q t |P S)2 ] (P s , 2 |<2t|P s ,i] 



(^ilQtl^] 



(C.3) 
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For n = 2 the result is 

ClQs, = ®! {^%^L(^ 3 + aAB + C) + {P a>1 - P s , 2 }} , (C.4) 

where we have defined 

a _ y ( 2 ^ - *0( 2i *J - gfj ~ 2 ^'(2Pj • K) A (2Q 3 -K)(2Q t ■ Q s ) - 2Q 2 (2Q t ■ K) 

{PsARj\Ps,2\ kPsAQt\Psl\ 

B = - k J2 ^ 2 Qs-K){2R r Q s )-2Q 2 s {2R J ■ K)f + 2Q^ 2 (p^p.jp^] (p s2 |p.|p sl ] 



j'=i 



(-Ps,l|-Rj|-fs,2] 



+ [(2Q S ■ XKgQt ■ Qs) - 2Q 2 (2Q t ■ K)} 2 - 2Q 2 K 2 (P a>1 \Qt\P a ,2] (P a ,2\Qt\P a ,i] 

c = ^ [(2Q a • Ar)(2^- • Q fl ) - 2gg (2^- • K)f - 3Q 2 K 2 (P^i^P^] (P s ,2\Rj\PsA 



J'=l 



(-Ps,l|-Rj|-Ps,2] 2 



2[(2Q S ■ XKggj ■ Q a ) - 2Q 2 (2R 3 ■ AT)] A 2[(2Q S ■ K)(2Q t ■ Qj - 2Q 2 (2Q t ■ K)] 

(PsMPsA {PsAQAPsA 

[(2Q S ■ K)(2Q t ■ Qs) - 2Q 2 (2Q t ■ K)} 2 - 3Q 2 S K 2 (P M |Q f |P a , 2 ] (P 8 , 2 \Qt\PsA 

(Ps,l\Qt\Ps,2\ 
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